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Qh Abstract 

We discuss some recent developments concerning the nucleon's helicity parton distribution func- 
, tions: New preliminary data from jet production at RHIC suggest for the first time a non- vanishing 

polarization of gluons in the nucleon. SIDIS measurements at COMPASS provide better constraints 
on the strange and light sea quark helicity distributions. Single-longitudinal spin asymmetries in 
^ 1^-boson production have been observed at RHIC and will ultimately give new insights into the 

light quark and anti-quark helicity structure of the nucleon. 

ON 
O 

2 1 Introduction 

^ QCD spin physics has been driven by the hugely successful experimental program of polarized deeply- 
inelastic lepton- nucleon scattering (DIS) \T\. One of the most important results has been the finding 
^ that the quark and anti-quark spins (summed over all flavors) provide only about a quarter of the 
nucleon's spin, AS ^ 0.25 in the proton helicity sum rule |2]: 



1 _ 1 

2 " 2 



= -AS + AG + L, + L, . (1) 



This result implies that sizable contributions to the nucleon spin should come from the gluon spin 
contribution AG, or from orbital angular momenta Lq^g of partons. To determine the other contributions 
to the nucleon spin has become a key focus of the field. In the present article, we describe some of 
the recent developments of the field. We focus on current efforts to determine the helicity parton 
distributions of the nucleon and on the latest experimental results. 

The helicity structure of the nucleon is foremost described by its twist-two helicity parton distribution 
functions, 

Af{x,Q^) = f+{x,Q^) - f-{x,Q^) {f = u,d,s,u,d,s,g), (2) 



*Talk presented by W. Vogelsang 



1 



/"•" (/~) denoting the number density of partons with the same (opposite) hehcity as the nucleon's, as a 
function of momentum fraction x and scale Q. QCD predicts the Q^-dependence of the densities through 
the spin-dependent Dokshitzer-Gribov-Lipatov-Altarelh-Parisi (DGLAP) evolution equations |3]: 



d f Aq 



dluQ^ \Ag 



\APgq{as,x) APgg{as,x) J \AgJ^ > 



(3) 



where ® denotes a convolution, and the sphtting functions AP^j are evaluated in QCD perturbation 
theory PSJE]. 

The contributions AS((5^) and AG((5^) in the helicity sum rule ([T]) are given by 



AS(g2) = / (At^ + Ai2 + A6i + A6i + As + As) (x, Q^)dx = AS^ + AS^ + AS, , 
Jo 



(4) 



AG{Q^) = I Ag{x,Q'')dx. (5) 







AS is independent of at the lowest order. The distributions have a proper field-theoretic definition. 
For example, in case of Ag it is given by [6] 

A^(x,g2) = _A_ I dXe^^^P^ {P,S\G+''{0)G\{Xn)\P,S)\^^ , (6) 

written in A'^ = gauge. Here, G^^ is the QCD field strength tensor, and G^^ its dual. The integral 
of Ag[x^ Q^) over all momentum fractions x becomes a local operator only in A'^ = gauge and then 
coincides with AG{Q^) [2j. 

The helicity parton distributions may be probed in spin asymmetries for reactions at large momen- 
tum transfer. The probes used so far are inclusive and semi-inclusive deep-inelastic lepton scattering 
(DIS and SIDIS, respectively), and pp scattering at large transverse momentum, see Fig. [l] Polarized 
DIS and SIDIS experiments have been carried out at SLAC, CERN, DESY and the Jefferson Laboratory 
P and mostly constrain the quark and anti-quark helicity distributions. RHIC at BNL [3, |8j is the 
first polarized proton-proton collider, operating at y/s = 200 and 500 GeV. The measurement of gluon 
polarization in the proton is a major focus and strength of RHIC. 

The basic theoretical concept that underlies much of spin physics is the factorization theorem. It 
states that large momentum-transfer reactions may be factorized into long-distance pieces that contain 




DIS SIDIS pp(RHIC) 

Figure 1: Parton- model Feynman diagrams for the processes constraining nucleon helicity structure. 
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the desired information on the spin structure of the nucleon in terms of its universal parton densities, 
and parts that are short-distance and describe the hard interactions of the partons. The latter can be 
evaluated using perturbation theory, thanks to the asymptotic freedom of QCD. As an example, we 
consider the reaction pp ttX, where the pion is produced at high transverse momentum p^, ensuring 
large momentum transfer. The statement of the factorization theorem [9j is then: 

d/\a = J] ® A/, ® dAal, ® (7) 

for the polarized cross section, where ® denotes a convolution. The D'l are the pion fragmentation 
functions. The sum in Eq. Q is over all contributing partonic channels a + 6 ^ c + X, with rfAa^^ 
the associated spin-dependent partonic cross section. Factorization is valid up to corrections that are 
suppressed as inverse powers of the hard scale. In general, a leading-order estimate of ([T]) merely captures 
the main features, but does not usually provide a quantitative understanding. Only with knowledge of 
the next-to-leading order (NLO) QCD corrections to the rfAa^^ can one reliably extract information on 
the parton distribution functions from the reaction. By now, NLO corrections are available for most of 
the processes relevant in polarized high-energy scattering [TO] . 

Independent information on the nucleon's helicity distributions may be obtained by using SU(2) and 
SU(3) flavor symmetries. The integrals of the flavor non-singlet combinations turn out to be proportional 
to the nucleon matrix elements of the quark non-singlet axial currents, (P, S \ qj^ 7^ g | P, S). Such 
currents typically occur in weak interactions, and one may relate the matrix elements to the /3-decay 
parameters F, D of the baryon octet. One flnds 

AS^ - ASrf = F + D = 1.267, 

AS^ + ASrf - 2AS, = 3F-D^ 0.58 . (8) 

If valid, the second relation when combined with Eq. Q gives that AS = 0.58 + SAS^, so that a small 
quark spin contribution to the proton spin implies a large negative strange quark contribution. Fairly 
signiflcant violations of SU(3) symmetry have been predicted based on heavy baryon chiral perturbation 
theory [llj. Lattice investigations of this issue have begun but are not yet conclusive [12j. 

2 Nucleon helicity structure: status 2009 

In recent publications [13j, we have presented the flrst next-to-leading order (NLO) "global" QCD 
analysis of the nucleon's helicity distribution from DIS, semi-inclusive DIS (SIDIS), and pp scattering 
at RHIC. We have used a Mellin moment method for the analysis. Our results are shown in Fig. [2} 
along with estimates of their uncertainties. The shaded bands in Fig. [2] show the distributions that are 
allowed if one permits an overall increase of Ax^ = 1 (green) or Ax^/x^ = 2% (yellow). As one can 
see, the "total" Au + Au and Arf + AJ helicity distributions are very well constrained. This is expected 
since these distributions are primarily determined by the large body of inclusive DIS data. Our results 
agree well with the distributions obtained in previous and other recent analyses [13 [lEl [T8] 
which considered only the lepton scattering data. 

The sea anti-quark distributions still carry rather large uncertainties, even though they are better 
constrained now than in previous analyses, thanks to the advent of more precise SIDIS data and of a 
new set of fragmentation functions [19] that describes the observables well in the unpolarized case. We 
flnd that the sea appears not to be SU(2)-flavor symmetric: the Au distribution is mainly positive, while 
the Ad anti-quarks carry opposite polarization. This pattern has been predicted at least qualitatively 
by a number of models [14^ 20j. Already based on the Pauli principle one would expect that if valence- 2i 
quarks primarily spin along the proton spin direction, uu pairs in the sea will tend to have the u quark 
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Figure 2: Present status of the nucleon's NLO helicity distributions according to the global analysis 
of Ref. [13]. The solid center lines show the best-fit result. The shaded bands provide uncertainty 
estimates, using a criterion of Ax^ = 1 (inner bands) or Ax^/x^ = 2% (outer bands) as allowed 
tolerance on the value of the fit. Also shown are results from earlier analyses [TUdS] of nucleon spin 
structure from lepton scattering data alone. 



polarized opposite to the proton. Hence, if such pairs are in a spin singlet, one expects Au > and, 
by the same reasoning. Ad < 0. We note that the uncertainties in SIDIS are still quite large, and it is 
in particular difficult to quantify the systematic uncertainty of the results related to the fragmentation 
mechanism at the relatively modest energies available. 

The strange sea quark density shows a sign change. At moderately large x ^ 0.1, it is constrained 
by the SIDIS data, which prefer a positive As. On the other hand, the inclusive DIS data combined 
with the constraints from baryon ;5-decays demand a negative integral of As. As a consequence. As 
obtains its negative integral purely from the contribution from low-x. Interestingly, there are initial 
lattice determinations of the integral AS^ [12], which point to small values. It is clearly important to 
understand the strange contribution to nucleon spin structure better. 

Constraints on the spin-dependent gluon distribution Ag predominantly come from RHIC. As can 
be seen from Fig. [2| the gluon distribution turns out to be small in the region of momentum fraction, 
0.05 < X < 0.2, accessible at RHIC, quite possibly having a node. At = 10 GeV^, the integral over 
the mostly probed x-region is found to be almost zero, /q 05 dxAg{x) = 0.005 ± 0.06, where the error is 
obtained for a variation of by one unit. Thus, on the basis of [13j, there are no indications of a sizable 
contribution of gluon spins to the proton spin. We also note that a way to access Ag in lepton-nucleon 
scattering at HERMES and COMPASS is to measure final states that select the photon-gluon fusion 
process, heavy-fiavor production and high-p^ hadron or hadron-pair production ^21] [22]. These data 
were not included in the analysis [13j, mostly because of the fact that success of the perturbative-QCD 
hard-scattering description had not been established for these observables in the kinematic regime of 
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Figure 3: COMPASS results [27] for SIDIS spin asymmetries on a proton target, compared to DSSV [13\ 
and DSSV+ fits [29j. 



interest here. At least for single-inclusive high-pr hadrons in h^X it has now been found that 

QCD hard-scattering does appear to be applicable in the COMPASS kinematic regime 



3 Recent developments 

Interesting new developments have taken place following the original DSSV analysis, mostly related to 
the advent of new data. We will summarize these in the following. 



3.1 Recent DIS and SIDIS data 

Recently, the COMPASS coUaboration has pubhshed new DIS [25] and SIDIS [26l [27] data. The latter 
extend the coverage in x down to about x 5 x 10~^, almost an order of magnitude lower than the 
kinematic reach of the HERMES data [28] used in the DSSV global analysis of 2008 [l3]. For the first 
time, the new results comprise measurements of identified pions and kaons taken with a longitudinally 
polarized proton target. Clearly, these data can have a significant impact on fits of helicity PDFs and 
estimates of their uncertainties. In particular, the new kaon data are expected to serve as an important 
check of the validity of the strangeness density obtained in the DSSV analysis discussed above, which 
instead of favoring a negative polarization as in most fits based exclusively on DIS data, prefers a 
vanishing or perhaps even slightly positive As in the measured range of x. 

Figure [3] shows a detailed comparison [29] between the new proton SIDIS spin asymmetries from 
COMPASS [27] and the original DSSV fit (dashed lines). Also shown is the result of a re- analysis 
at NLO [29] accuracy (denoted as "DSSV+") based on the updated data set. The differences between 
the DSSV and the DSSV+ fits are hard to notice, both for identified pions and kaons. The total 
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Figure 4: profiles for the truncated first moment of As in two different x intervals, 0.02 < x < 1 
(left) and 0.001 <x < 0.02 (right). 



of the fit drops only by a few units upon refitting, which is not really a significant improvement for a 
PDF analysis in view of non-Gaussian theoretical uncertainties. The change in is also well within 
the maximum Ax^/x^ = 2% criterion adopted in the original DSSV global analysis [13]. Overall, upon 
refitting, one finds a trend towards smaller net polarization for Au and Ad in the range 0.001 < x < 1 
than in DSSV. 

As we saw earlier, the original DSSV fit fT3] found an interesting feature for the strangeness helicity 
distribution: As was found to be small and slightly positive at medium-to-large x, but has a significantly 
negative first moment in accordance with expectations based on SU(3) symmetry and fits to DIS data 
only. To investigate this issue further, we present in Fig. [4] the x^ profiles for two diflterent intervals 
in x: 0.02 < x < 1 (left) and 0.001 < x < 0.02 (right). The profiles in Fig. [i] clearly show that the 
result for As for 0.001 < x < 0.02 is a compromise between DIS and SIDIS data, the latter favoring 
less negative values. Interestingly though, the new COMPASS SIDIS data, which extend towards the 
smallest x values so far, actually show some preference for a slightly negative value for As as well. 
For 0.02 < X < 1 everything is determined by SIDIS data, and all sets consistently ask for a small, 
slightly positive strange quark polarization. There is no hint of a tension with DIS data here as they 
do not provide a useful constraint at medium-to-large x. We note that at low x most SIDIS sets in the 
original DSSV fit give indifferent results. We also mention that in the range x > 0.001 the hyperon 
decay constants, the so-called F and D values, do not play a significant role in constraining As{x). 
To quantify possible SU(3) breaking effects one needs to probe As{x) at even smaller values of x, for 
instance in SIDIS at a future EIC t3Qj. We finally note that the HERMES and COMPASS data are 
consistent in the region of overlap, 0.02 < x < 1. 

Clearly, all current extractions of As from SIDIS data suffer from a significant dependence on kaon 
FFs, see, e.g., Refs. lESJEZ], and better determinations of D^[z) are highly desirable. Contrary to other 
fits of FFs [32j, only the DSS sets [19j provide a satisfactory description of pion and kaon multiplicities 
in the same kinematic range where we have polarized SIDIS data. 
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3.2 W bosons at RHIC 



We have seen in the previous section that the SIDIS data provide some insights into the flavor structure 
of the polarized sea distributions of the nucleon, albeit with still fairly large uncertainties. Complemen- 
tary and clean information on At^, AtZ, Arf, and Arf will come from pp W^X at RHIC, where one 
exploits the maximally parity- violating couplings of the produced W bosons to left-handed quarks and 
right-handed anti-quarks [71 [33] . These give rise to a single-longitudinal spin asymmetry, 

for the processes pp £^X^ where the arrow denotes a longitudinally polarized proton and ^ = e or 
fi is the charged decay lepton. The high scale set by the W boson mass makes it possible to extract 
quark and anti-quark polarizations from inclusive lepton single-spin asymmetries in W boson production 
with minimal theoretical uncertainties, as higher order and sub-leading terms in the perturbative QCD 
expansion are suppressed [3H |35l [36l |37j . 

For W~ production, neglecting all partonic processes but the dominant ud W~ one, the spin- 
dependent cross section in the numerator of the asymmetry is found to be proportional to the combi- 
nation 

Act oc Au{xi) (8) rf(x2)(l — cos^)^ — Ad{xi )®^(x2)(l + cose)^ (10) 

where 9 is the polar angle of the electron in the partonic c.m.s., with 9 = in the forward direction of 
the polarized parton. At large negative pseudorapidity ryiept of the charged lepton, one has X2 ^ Xi and 



9 ^ 7t/2. In this case, the first term in Eq. (10) strongly dominates, since the combination of parton 
distributions, A'u(xi)rf(x2), and the angular factor, (1 — cos^)^, each dominate over their counterpart 
in the second term. Since the denominator of is proportional to u{xi) rf(x2)(l — cos^)^ + d{xi) 
u(x2){l + cos^)^, the asymmetry provides a clean probe of Au(xi)/u{xi) at medium values of Xi. By 



similar reasoning, at forward rapidity r^iept ^ the second term in Eq. (10) dominates, giving access to 
— Arf(xi)/rf(xi) at relatively high xi. 

For production, one has the following structure of the spin-dependent cross section: 

Act oc Ad{xi) u{x2){l + cos9f - Au{xi) ® d{x2){l - cos9f . (11) 

Here the distinction of the two contributions by considering large negative or positive lepton rapidities is 
less clear-cut than in the case of W~ . For example, at negative ryiept the partonic combination d{xi)u{x2) 
will dominate, but at the same time 9 ^ 7r/2 so that the angular factor (1 + cos 6^)^ is small. Likewise, 
at positive ryiept the dominant partonic combination Au{xi)d{x2) is suppressed by the angular factor. 



So both terms in Eq. (11) will compete essentially for all r^iept of interest. Our global analysis technique 
is of course suited for extracting information on the polarized PDFs even if there is no single dominant 
partonic subprocess. The NLO corrections to the single-inclusive lepton cross sections have recently 
been presented in [37] in a way tailored to use in the global analysis framework. 

Figure [5] shows the first published data from RHIC for Al in production [MlESj- For now, the 
statistical uncertainties are still large. However, already now a large negative asymmetry is seen for the 
case of production, resulting primarily from the positive up-quark polarization in the proton (see 
Eq. ([TT])). Clearly, there is a large potential in future Vl^-measurements at RHIC. 

3.3 New constraints on Ag 

The STAR and PHENIX experiments at RHIC have recently presented new preliminary data from 
the 2009 run for the double-helicity spin asymmetry All for jet and neutral-pion production, respec- 
tively ^401141]. The results are shown in Figs. |6] and [Tj One can see that the experimental uncertainties 
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Figure 5: Published STAR [38j (top) and PHENIX [39j (bottom) data for the single-hehcity asymmetry 
Aliu production at RHIC. 



are very significantly reduced as compared to those in the previous run-6 data sets [ill |43]. An inter- 
esting feature of the new preliminary STAR data is that they lie consistently above the result for the 
best-fit DSSV distribution for jet transverse momenta below 25 GeV or so. They do remain well below 
the old GRSV- "standard" result of fT4], on the other hand. This suggests that the spin-dependent 
gluon distribution may be somewhat different from zero in the x-range where it is constrained by the 
RHIC data. The thick (magenta) solid line in the figure shows All obtained for a special set of par- 
ton distributions within the DSSV analysis. For this set the truncated moment of /\g over the region 
0.001 < X < 1 was varied, allowing the total to change by 2%. Evidently, this set of parton distribu- 
tions describes the STAR data rather well. The truncated moment of /S.g in this set over the x-range 
accessed at RHIC is 

/ dx Ag{x, = 10 GeV^) = 0.13 , (12) 
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which is just within the range 

rO.2 

/ dx Ag{x, = 10 GeV^) = 0.0051^;}^^ (13) 

^0.05 

quoted as more conservative uncertainty {Ax^/x^ = 2%) in [13\. Despite the fact that this really is 
only an illustration that cannot replace a proper re-analysis of the data, it does appears that, for the 
first time, there are indications of non- vanishing gluon polarization in the nucleon. Figure [7| shows the 
comparison to the new preliminary PHENIX tt^ data Here the values of All are much smaller, 
which is mostly due to the fact that lower values of x are probed at the transverse momenta relevant in 
the PHENIX measurements. One can see that the data are well described by both the DSSV set and 
the special set of polarized parton distributions used in Fig. |6| 

4 The future: Electron Ion Collider (EIC) 

An Electron Ion Collider is currently being considered in the U.S. as a new frontier facility to explore 
strong-interaction phenomena [31j. One of its key features would be the availability of high-energy, 
high-luminosity polarized ep collisions to probe nucleon spin structure. This would also allow precision 
extractions of Ag^ in particular from scaling violations of the proton's spin-dependent structure function 
gi. Figure [S] shows the results of a recent dedicated phenomenological study [44j. "Pseudo" ElC-data 
were generated for collisions of 5 GeV electrons with 50, 100, 250, and 325 GeV protons and were added 
to the DSSV global analysis code. The statistical precision of the data sets for 100 — 325 GeV protons 
was taken to correspond to about two months of running at the anticipated luminosities for eRHIC 
with an assumed operations efficiency of 50%. For 5 x 50 GeV an integrated luminosity of 5fb~^ was 
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Figure 7: As in Fig. pi but for the preliminary PHENIX run-9 data ^T] . 



assumed. The projected uncertainties were used to randomize the pseudo-data by one sigma around 
their central values determined by the DSSV set of PDFs. With these projected EIC data with their 
estimated uncertainties, a re-fit of the DSSV polarized parton distribution functions was performed. 
The results are shown in Fig. [8} As one can see, with EIC data it should be possible to map the 
currently completely undetermined shape of Ag for 10"*^ ^ x < 0.01 to an accuracy of about ±10% or 
better. 
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